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ABSTRACT 

We use the ages, masses and metallicities of the rich young star cluster systems in the 
nearby starburst galaxies NGC 3310 and NGC 6745 to derive their cluster formation 
histories and subsequent evolution. We further expand our analysis of the systematic 
. uncertainties involved in the use of broad-band observations to derive these param- 

■^j- ' eters (Paper I) by examining the effects of a priori assumptions on the individual 

OO ! cluster metallicities. The age (and metallicity) distributions of both the clusters in the 

circumnuclear ring in NGC 3310 and of those outside the ring are statistically indis- 
lO ' tinguishable, but there is a clear and significant excess of higher-mass clusters in the 

, ring compared to the non-ring cluster sample; it is likely that the physical conditions 

in the starburst ring may be conducive for the formation of higher-mass star clusters, 
on average, than in the relatively more quiescent environment of the main galactic 
r"| . disc. For the NGC 6745 cluster system we derive a median age of ~ 10 Myr. NGC 

6745 contains a significant population of high-mass "super star clusters" , with masses 
.1 | in the range 6.5 < log(M c i/Af ) < 8.0. This detection supports the scenario that 

such objects form preferentially in the extreme environments of interacting galaxies. 
The age of the cluster populations in both NGC 3310 and NGC 6745 is significantly 
, lower than their respective characteristic cluster disruption time-scales, respectively 

log(i4 ls /yr) = 8.05 and 7.75, for 10 4 M Q clusters. This allows us to obtain an indepen- 
dent estimate of the initial cluster mass function slope, a = 2.04(±0.23)±^ for NGC 
3310, and 1.96(±0.15) ± 0.19 for NGC 6745, respectively, for masses M d > 10 5 M Q 
and M c \ >4x 10 5 M Q . These mass function slopes are consistent with those of other 
young star cluster systems in interacting and starburst galaxies. 

Key words: galaxies: individual: NGC 3310, NGC 6745 - galaxies: starburst - 
galaxies: star clusters 



1 INTRODUCTION 

1.1 Star cluster formation in intense starbursts 

The production of luminous, compact star clusters seems to 
be a hallmark of intense star formation. Such young clusters 
have been identified in intense starburst regions in several 
dozen galaxies, often involved in interactions (e.g., Holtz- 
man et al. 1992, Whitmore et al. 1993, O'Connell et al. 
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1994, Conti et al. 1996, Watson et al. 1996, Carlson et al. 
1998, de Grijs et al. 2001, 2003d). Their sizes, luminosities, 
and - in several cases - spectroscopic mass estimates are 
entirely consistent with what is expected for young Galaxy- 
type globular clusters (GCs; Meurer 1995, van den Bergh 

1995, Ho & Filippenko 1996a,b, Schweizer & Seitzer 1998, 
de Grijs et al. 2001, 2003a). 

It is possible, even likely, that a large fraction of the 
star formation in starbursts takes place in the form of such 
concentrated clusters (see, e.g., de Grijs et al. 2001, 2003d 
for a discussion). Young compact star clusters are there- 
fore important because of what they can tell us about GC 
formation and evolution (e.g., destruction mechanisms and 
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efficiencies). They are aiso important as probes of the his- 
tory of star formation, chemical evolution, initial mass func- 
tion (IMF), and other physical characteristics in starbursts. 
This is so because each cluster approximates a coeval, single- 
metallicity, simple stellar population (SSP). Such systems 
are the simplest to model, after individual stars themselves, 
and their ages and metallicities and - in some cases - IMFs 
can be estimated from their integrated spectra. 

In de Grijs et al. (2003c; hereafter Paper I) we developed 
a reliable method to determine cluster ages, masses, metal- 
licities and extinction values robustly and simultaneously 
based on imaging observations in a minimum of four broad- 
band passbands covering the entire optical wavelength range 
from the U to the / band or their equivalents in non-standard 
passband systems. We tested our method using the ~ 150 
clusters in the centre of the nearby starburst galaxy NGC 
3310, and confirmed the previously suggested scenario that 
NGC 3310 underwent a (possibly extended) global burst of 
cluster formation ~ 3 x 10 7 yr ago, likely triggered by the 
last tidal interaction or merger with a low-metallicity, gas- 
rich dwarf galaxy (see Paper I). 

However, in determining the cluster formation history 
from the age distribution of magnitude-limited cluster sam- 
ples, cluster disruption must be taken into account. This is 
because the observed age distribution is that of the surviv- 
ing clusters only. In this paper, we will therefore embark on 
a further analysis of the formation history of the NGC 3310 
cluster system and its subsequent evolution. We will also 
apply our method to the star clusters and star-forming re- 
gions in the interacting galaxy NGC 6745 (the "Bird's Head 
Galaxy"), and interpret its cluster formation history in the 
contexts of its recent tidal encounter and of cluster disrup- 
tion. 

1.2 An empirical description of cluster disruption 

As shown by Boutloukos & Lamers (2003, hereafter BL03), 
and applied to the fossil starburst cluster sample in M82 
B by de Grijs et al. (2003b), with only a few well-justified 
assumptions, the mass and age distributions of a magnitude- 
limited sample of star clusters in a given galaxy can be pre- 
dicted both accurately and robustly, despite the complex 
physical processes underlying the assumptions (for a full dis- 
cussion see BL03). If all of the following conditions are met 
then it can be shown easily that the age distribution of the 
observed cluster population will obey approximate power- 
law behaviours. 

(i) The cluster formation rate, dN d f c1 ' 1 — S x M~^ a , is 
constant; 

(ii) The slope a of the cluster IMF is constant among 
cluster systems, with N{M C \)&M oc M~ x a AM; 

(iii) Stellar evolution causes SSPs to fade as F\ ~ t^^ , as 
predicted by theoretical cluster evolution models; and 

(iv) The cluster disruption time depends on their initial 
mass as t dis = tf s x (M c i/10 4 M G ) 7 , where tf B is the disrup- 
tion time-scale of a cluster with initial mass M c i = 1O 4 M . 

It is well-established, however, that the disruption time-scale 
does not only depend on mass, but also on the initial cluster 
density and internal velocity dispersion (e.g., Spitzer 1957, 
Chernoff & Weinberg 1990, de la Fuente Marcos 1997, Porte- 
gies Zwart et al. 2001). Following the approach adopted in 



BL03 and de Grijs et al. (2003b), we point out that if clus- 
ters are approximately in pressure equilibrium with their 
environment, we can expect the density of all clusters in a 
limited volume of a galaxy to be roughly similar, so that 
their disruption time-scale will predominantly depend on 
their (initial) mass (with the exception of clusters on highly- 
eccentric orbits) . In the opposite case that the initial cluster 
density p depends on their mass M c i in a power-law fashion, 
e.g., p ~ Af°j with x being the (arbitrary) power-law expo- 
nent, the disruption time-scale will also depend on mass if 
t dis ~ M C V (BL03). 

The observed cluster age distribution will obey the fol- 
lowing approximate power-law behaviours: 

• diVd/di oc t'>' 1 ~ Q1 " ) for young clusters as a result of fad- 
ing; 

• diVd/di oc f t 1- ™)/""' for old clusters as a result of disrup- 
tion. 

Similarly, the mass spectrum of the observed clusters will be 

• diVci/dMd oc M C ( 1 1/C)_Q for low- mass clusters as a result 
of fading; 

• diV c i/dM c i oc M^~ a for high- mass clusters as a result 
of disruption. 

Thus, both distributions will show a double power law 
with slopes determined by a, ( and 7. The crossing points, 
t cr oss and M cross , are determined by the cluster formation 
rate and by the £4" time-scale. In these expressions, N c i 
is the total number of clusters in a given sample, and so 
diV c i/dM c i and diV c i/d£ are the numbers of clusters per con- 
stant bin in mass and age, respectively. 

In BL03 and de Grijs et al. (2003b) we showed that 
the observed age and mass distributions of the star clus- 
ter systems in a number of well-studied galaxies indeed 
show the predicted double power-law behaviour, where ap- 
plicable modified by a non-constant (bursty) cluster forma- 
tion rate. From the analysis of these observed distributions 
BL03 showed that the value of 7 is approximately constant, 
(7 = 0.62 ± 0.06), under the very different environmen- 
tal conditions in their sample galaxies, but the characteris- 
tic disruption time-scales (i.e., t^) differ significantly from 
galaxy to galaxy. 

In Sect. 2 we will derive the formation history of the 
NGC 3310 cluster system; we will distinguish between star 
clusters in the starburst ring and those detected outside 
(Sect. 2.1), and interpret our results in the context of star 
cluster disruption processes and time-scales (Sect. 2.2). We 
will then discuss the system of young star clusters and star- 
forming regions in NGC 6745 in Sect. 3 in the context of its 
very recent tidal encounter with a small companion galaxy 
(Sect. 3.2), and explore the cluster properties in the frame- 
work of their disruption history (Sections 3.3 and 3.4). As we 
will see, both cluster systems are too young to have already 
undergone significant disruption, so that we are in fact ob- 
serving their properties in close to initial conditions. We will 
discuss the derived slopes for the initial cluster mass func- 
tions for both galaxies in the general context of star cluster 
formation in Sect. 4, and then summarize our results and 
conclusions in Sect. 5. 
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2 MEASURING THE INITIAL CLUSTER 
MASS FUNCTION SLOPE IN NGC 3310 

In Paper I we obtained robust age and mass estimates for 
~ 150 star clusters in the nearby starburst galaxy NGC 
3310. We used a number of passband combinations based on 
archival Hubble Space Telescope (HST) observations from 
the ultraviolet (UV) to the near-infrared (NIR) to achieve 
this. In this paper, we will use these age and mass estimates 
to derive further details of the clusters' formation history 
and their subsequent evolution. 

Figure la shows the distribution of the NGC 3310 clus- 
ters in the age vs. mass plane. The solid line overplotted on 
the figure shows the expected effect of evolutionary fading 
of an instantaneously formed SSP. The fading line shown is 
based on zero extinction, and applies to an observed cluster 
system with a limiting magnitude of V = 21 at the dis- 
tance of NGC 3310 (m — M = 30.57; see Paper I), and our 
adopted, Salpeter-type IMF (with stellar masses m in the 
range 0.15 < tu/Mq < (50 — 70), the upper limit depend- 
ing on the metallicity and determined by the mass coverage 
of the Padova isochrones) , predicted by the Gottingen SSP 
models (Schulz et al. 2002, updated to include nebular emis- 
sion by Anders & Fritze-v. Alvensleben 2003; see also Paper 
I). The predicted lower limit agrees very well with our data 
points. 

Because of the way in which we performed our cluster 
selection in Paper I, our adopted photometric completeness 
limit is ultimately determined by the interplay of two effects. 
These include the photometric depth of the shallowest ex- 
posure for which we required genuine source detections (i.e., 
the F300W image, for which we determined an average 90% 
completeness fraction at V — 22; sec Figure 1 in Paper I) 
and the further loss of completeness due to the final, visual 
verification of the cluster candidates detected automatically. 
We conservatively estimate this latter step, which dominates 
the V-band photometry, to reduce our sample completeness 
by of order one more magnitude (see the detailed discussion 
on completeness effects for the analysis of the M51 cluster 
system in N. Bastian et al., in prep.), so that the limiting 
magnitude at ~ 90% completeness for the full sample of 
NGC 3310 clusters is V ~ 21. 

For a nominal extinction of Ay = 0.2 mag (see Pa- 
per I), the detection limit will shift to higher masses by 
Alog(M c i/Mo) ~ 0.08, which is well within the uncertain- 
ties associated with our mass determinations. 

Figures lb and c show the distributions of, respectively, 
the ages and masses of the compact clusters in the centre 
of NGC 3310. The uncertainties in our fit results are of the 
order of the histogram bin sizes (Paper I). The open his- 
tograms, with their associated Poissonian error bars, repre- 
sent the full cluster sample, whereas the shaded histograms 
show the effects of adopting a fixed age or mass cut-off. 
The shaded histogram in Fig. lb corresponds to the clusters 
with masses greater than log(M c i/MQ) = 5.0, as indicated 
by the horizontal dotted line in Fig. la (arbitrarily chosen), 
and the shaded histogram in Fig. lc represents (for an arbi- 
trary age cut-off) the youngest clusters, log(Age/yr) < 7.65, 
corresponding to the data points to the left of the vertical 
dotted line in Fig. lc. The age and mass-dependent com- 
pleteness limit introduces some skewness into the age and 
mass distributions compared to fixed age and mass cut-offs. 



2.1 Starburst ring versus non-ring clusters 

The morphology over a wide range of wavelengths (from 
X-rays to radio waves) of the central regions of NGC 3310 
is dominated by a bright, relatively dense ring- like struc- 
ture containing a large number of actively star-forming re- 
gions and young star cluster candidates (see Paper I for a 
review). Conselice et al. (2000) suggested a bar-driven ori- 
gin for the active starburst in the ring, combined with the 
recent infall of a companion galaxy. Such a scenario provides 
a natural explanation for the low metallicity observed in the 
star-forming knots near the galactic centre (see Paper I for 
an overview), while it also explains why we observe concen- 
trated star formation in star clusters or luminous Hn regions 
in such a tightly-wound ring-like structure surrounding the 
centre (e.g., van der Kruit & de Bruyn 1976, Telesco & Gat- 
ley 1984, Pastoriza et al. 1993, Meurer et al. 1995, Smith 
et al. 1996, Conselice et al. 2000, Elmegreen et al. 2002), 
coinciding with the end of the nuclear bar (Conselice et al. 
2000, but see Diaz et al. 2000), but not inside this ring. 

2.1.1 A comparison of basic properties 

We identified 82 of the 147 star clusters in our sample to 
coincide with this circumnuclear ring, with the 65 remain- 
ing objects being located outside this starburst ring. In Fig. 
2 we compare the basic properties of these cluster sam- 
ples. The peaks in the age distributions are observed at 
log(Age/yr) = 7.45 and 7.43 for the clusters in and out- 
side the ring, respectively, with corresponding Gaussian sig- 
mas, o~ g , of 0.40 and 0.38 in logarithmic age space. How- 
ever, while the age (and also the metallicity) distributions 
of both cluster samples are statistically indistinguishable, 
the ring clusters appear to peak at slightly greater masses 
than those found outside the ring: the mass distribution of 
the clusters in the ring peaks at log(M c i/Mo) r i ng = 5.29, 
versus log(M c i/MQ) n0 n-ring = 5.13 for the clusters outside 
the ring. Their Gaussian sigmas are also significantly differ- 
ent, at <7G,rmg = 0.32 and ctg, non-ring = 0.51, respectively. 

We examined whether this difference could be due to a 
spatial dependence of the completeness fraction, in the sense 
that one might expect a lower completeness in the ring area 
with its higher background flux than outside the ring. This 
difference in the mean (logarithmic) mass corresponds to 
a difference in the peak of the cluster luminosity function 
(CLF), in any passband, of Am pe ak — 0.4 mag, and the dif- 
ference in the Gaussian sigmas of the distributions is also 
in the sense expected if this were due to different and vari- 
able levels of completeness. We examined the mean surface 
brightness levels in and outside the starburst ring, respec- 
tively, both in the F606W image and in the slightly shallower 
F300W passband. The difference in the mean surface bright- 
ness levels between the starburst ring and the area outside 
the ring is A(£ n0 n-ring — £rm g ) ~ 0.3 — 0.9 mag arcsec -2 
and ~ 0.1 - 0.8 mag arcsec~ 2 in the F300W and F606W 
images, respectively, based on a per-pixel comparison. 

The overabundance of low-mass clusters outside the 
ring compared to their counterparts in the starburst ring 
might therefore (at least partially) be due to a spatial de- 
pendence of the completeness fraction. However, there is a 
clear and significant excess of higher-mass clusters (5.25 < 
log(M c i/M ) < 5.75) in the ring compared to the non-ring 
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Figure 1. (a) - Distribution of the NGC 3310 clusters in the (age vs. mass) plane. Overplotted is the expected, age-dependent detection 
limit. This model prediction is based on a detection limit of rrtpgoew ~ V = 21 (see text) and (m — M)ngc3310 = 30.57, assuming no 
extinction. The features around 10 Myr are caused by the appearance of red supergiants. (b) and (c) - Distributions of, respectively, the 
ages and masses of the compact clusters in the centre of NGC 3310. The open histograms, with their associated Poissonian error bars, 
represent the full cluster sample; the shaded histogram in panel (b) corresponds to the clusters with masses greater than log(Af c i/Afo) = 
5.0, indicated by the horizontal dotted line in panel (a), and the shaded histogram in panel (c) represents the youngest clusters, 
log(Agc/yr) < 7.65, corresponding to the data points to the left of the vertical dotted line in panel (a). The uncertainties arc of the order 
of the histogram bin sizes. 



cluster sample, even after taking the systematic uncertain- 
ties in our mass estimates into account (the estimated ~ la 
systematic uncertainties are of the order of the histogram bin 
sizes). If both the ring and the non-ring cluster populations 
were drawn from the same parent population, one would not 
expect to observe such large differences, and thus it appears 
that the physical conditions in the starburst ring, such as 
caused by the effects of the proposed bar-driven instabili- 
ties, the higher density of the interstellar medium (ISM), 
and the associated higher likelihood to encounter significant 
propagating shock waves, may be conducive for the forma- 
tion of higher-mass star clusters, on average, than in the 
relatively more quiescent environment of the main galactic 
disc. Thus, the origin of the differences in the global mass 



distributions at the high-mass end of the two cluster samples 
is most likely found in their respective ISM properties. 

If the cluster formation in the starburst ring were pre- 
dominantly due to the bar-driven instabilities suggested by 
Conselice et al. (2000), one would, to first order, expect a 
narrower age distribution for the clusters in the ring with 
respect to that of the non-ring clusters in the general field 
of the central galactic disc. However, both age distributions 
are characterized by a median age of (log(Age/yr)) ~ 7.5 
(~ 30 Myr), and an age spread of aG(log[Age/yr]) ~ 0.4 
(~ 70 Myr, from ~ 80 — 10 Myr ago). Furthermore, if we 
examine the age distributions of the clusters in linear age 
space, it appears that clusters both in and outside the star- 
burst ring have been forming approximately constantly, but 



© 2003 RAS, MNRAS 000, 000-000 



Star Cluster Formation and Evolution in Nearby Starburst Galaxies: II. 5 



NGC 3310 




6 7 8 9 10 
log( Age / yr ) 




4 5 6 7 
log( M ci / M ) 




0.02 0.04 0.06 
Metallicity Z 



E^^sl Clusters in ring (82) V////A Clusters outside ring (65) 



Figure 2. Comparison of the basic properties of the ring clusters with those of the non-ring cluster sample. 



at a ~ 2 — 3x higher level during the past ~ 40 Myr than 
before. Thus, we conclude that the starburst producing the 
NGC 3310 ring clusters has been ongoing for at least the 
past 40 Myr, at an approximately constant cluster forma- 
tion rate during the burst period. Star cluster formation 
has proceeded at a similar rate in the general central disc 
of the galaxy; due to the extreme youth of both samples of 
starburst-induced star clusters, statistically significant dif- 
ferences between the age distributions of the ring and non- 
ring populations have not yet had time to develop, which 
implies that we are still observing at least part of the clus- 
ter sample in the environment determined by their initial 
conditions. 



2.1.2 Effects of metallicity changes on the systematic 
uncertainties 

Finally, we caution that the peak of cluster formation, while 
relatively robust (see, e.g., Paper I for a discussion), is un- 
doubtedly broadened by the various uncertainties entering 
the age-dating process. Some of these (systematic) uncer- 
tainties were discussed in detail in Paper I, where we con- 
centrated our discussion on the choice of passband combi- 
nation for this type of analysis. Here, we will try to obtain a 
further handle on the systematic uncertainties. If the cluster 
formation in the starburst ring were induced by bar-driven 
instabilities acting on very short time-scales (i.e., on time- 
scales short compared to the median age of the cluster sam- 
ple), one would not expect the resulting cluster sample to be 
characterized by a large range in metallicity. This argument 
holds both if the progenitor gas clouds originated from the 
main NGC 3310 gas reservoir, and also if they came from a 
low-metallicity dwarf galaxy that was cannibalized (assum- 
ing that such a galaxy would exhibit a fairly small range 
in metallicities, as suggested for this particular interaction; 
see Paper I). The former case is further strengthened by ob- 
servations in other galaxies that old galactic discs generally 
show smooth, very slowly declining radial metallicity gradi- 
ents, such that at any given small distance range from the 
galactic centre (e.g., as in the NGC 3310 starburst ring) the 
stellar metallicity is approximately constant (for the Galaxy 



see, e.g., Twarog et al. 1997, and the discussion therein). 
This is to some extent already reflected in Fig. 2c, where we 
derived that > 80% of the clusters in the ring have metal- 
licities Z < 0.015. 

We will now assume a mean metallicity of Z — 0.008 
for the clusters in the starburst ring, and rederive their age 
and mass distributions. This will give us a further handle 
on the systematic, model-dependent uncertainties entering 
the age-dating and mass-estimation process. In Fig. 3 we 
show the age and mass distributions obtained under the as- 
sumption that the ring clusters all have a fixed metallicity 
of Z = 0.008. The peaks in the age and mass distributions 
occur at (log(Age/yr)) = 7.34 and (log(M c i/M Q )) = 5.09, 
respectively, while the distributions are characterized by 
Gaussian sigmas o"G.ag C = 0.50 and o"G,mass = 0.42. As a 
reminder, the age and mass distributions resulting from re- 
taining the metallicity as a free parameter (Sect. 2.1.1), were 
(log(Age/yr)> ~ 7.45 and (log(M c i/M Q )) = 5.29, respec- 
tively, with <7G,a g c = 0.40 and (TG.mass = 0.32. We see that 
while the peak values are retained with a reasonable ro- 
bustness (although there is, as expected, some effect caused 
by the age—met allicity degeneracy), both distributions have 
broadened significantly. This confirms, therefore, that (i) it 
is very important to determine all of the free parameters 
(age, metallicity, and extinction, and the corresponding mass 
estimates) independently for each individual cluster, instead 
of assuming a generic value for any of these parameters (if a 
sufficient number of broad-band passbands covering a large 
wavelength range are available); and (ii) the widths of the 
distributions (assuming a fixed metallicity) are broader than 
their intrinsic widths owing to the propagation of model and 
measurement uncertainties (e.g., de Grijs et al. 2003b). 

For comparison, we also computed the characteristics 
of the age and mass distributions for the NGC 3310 clus- 
ters outside the ring, under the same assumption of a fixed 
metallicity of Z — 0.008 for all clusters. The trend between 
these new determinations and those of the fixed-metallicity 
results for the ring clusters remains similar as for the situa- 
tion in which we allowed the metallicity to be an additional 
free parameter: the age distribution peaks again at a similar 
median age, (log(Age/yr)) non _ r i ng = 7.37, with a Gaussian 
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Figure 3. Age and mass distributions of the clusters in the NGC 
3310 starburst ring, obtained under the assumption that all clus- 
ters have a fixed mctallicity of Z = 0.008. 



sigma of a G ,, 



= 0.52. The median of the mass dis- 



tribution occurs at a significantly lower mass than for the 
ring clusters, (log(M c i/MQ)} n0 n-ring = 4.93, with a Gaus- 
sian sigma of a G , non-ring = 0.60. Thus, we arrive at similar 
conclusions as before for the mass distribution of the non- 
ring compared to the ring sample of clusters in NGC 3310. 



2.2 The star cluster disruption context 

In Fig. 4 we show the formation rate and the mass spectrum 
of the entire NGC 3310 cluster sample. These distributions 
depend on the cluster formation history and on the cluster 
disruption time-scale governing the centre of NGC 3310. 

For our analysis of the cluster disruption history in the 
centre of NGC 3310 we adopt a slope C = 0.648 for the 
evolutionary fading of clusters in the V band, derived from 
theoretical SSP models (see also BL03, de Grijs et al. 2003b). 
For the mass scaling of the disruption time-scale we adopt 
7 = 0.62 (BL03, de Grijs et al. 2003b). 

We also assume that the cluster formation rate has been 
approximately constant during the periods of interest for 
our analysis, log(Age/yr) < 7.5 and log(Age/yr) > 7.5. As 
shown in Sect. 2.1, this assumption seems justified, within 
the uncertainties. 

In Fig. 4a, the long-dashed line for the youngest ages, 
log(Age/yr) < 8, is the expected effect of a fading SSP 
(in the V band), for a constant ongoing cluster formation 
rate, shifted vertically to best match the data points for 
log(Agc/yr) < 7.5. The short-dashed lines in both Figs. 4a 
and b are the best-fitting power-law slopes, obtained for the 
age and mass ranges log(Age/yr) > 7.7 and log(M c i/M ) > 
5.1, while the dotted lines are indicative of the likely un- 
certainties in the slopes. The latter were obtained from 
subsets of the data points used for the dashed fits, deter- 
mined over the full ranges. Our results are robust with re- 
spect to changes in the adopted bin size in log(Age/yr) 
and log(M c i/MQ). The vertical dashed line in Fig. 4b at 
\og(M c \/ Mq) — 4.4 corresponds to our completeness limit, 
for a cluster age of log(Age/yr) = 7.5. However, since our 
cluster sample is characterized by a range of ages up to ~ 10 8 
yr, it follows that incompleteness plays at least some role for 
masses up to log(M c i/M Q ) ~ 4.9 (see the vertical dotted line 
in Fig. 4b). For the age distribution the effects of incomplete- 
ness are less complicated, since they are determined solely 
by (the extension of) the fading line in Fig. 4a. Moreover, 
the effects of the non-constancy of the cluster formation rate 



are less transparent and therefore harder to disentangle from 
the mass distribution than from the age distribution (see de 
Grijs et al. 2003b for a discussion). 

The short-dashed line in Fig. 4a, which shows the effect 
of cluster disruption, has a predicted slope of (1 — a)/7 and 
an observed slope of — 2.64(±0.14) ± 0.47, where the first 
uncertainty represents the formal uncertainty in the fit and 
the second the likely uncertainty resulting from our choice 
of fitting range (i.e., the dotted fits in Figs. 4a and b). If we 
adopt a cluster IMF slope of a — 2.0 (see Sect. 4) we find 
that the slope of the cluster disruption law is 7 ~ 0.38 — 0.46, 
with £ d is oc MJj. This value is smaller than the mean value 
of 7 ~ 0.62 found by BL03 from cluster samples in differ- 
ent galaxies. The crossing point of the two power-law fits 
is at logt cr oss — 7.7. Substituting this value in Eq. (15) of 
BL03, we find that the characteristic disruption time of a 
1O 4 M cluster is logtf 3 = 7.9. This implies that the dis- 
ruption time-scale of clusters with an initial mass greater 
than 1O 5 M will be greater than 2 x 10 8 yr. The slope of 
the mass distribution in Fig. 4b, for clusters more massive 
than 10 5 M Q is -2.04(±0.23)tg where the first uncer- 
tainty represents the formal uncertainty in the fit and the 
second the likely uncertainty resulting from our choice of fit- 
ting range. All these clusters are younger than 2 x 10 yr. (In 
fact, all our sample clusters, except one, are younger than 
this age.) Thus, disruption has not yet affected the mass dis- 
tribution of these massive clusters and hence the observed 
slope of the mass distribution is the slope of the cluster IMF. 
We note that this new value of the NGC 3310 cluster IMF 
slope closely matches both our earlier CLF slope determina- 
tion of a — 1.8 ±0.4 obtained from the cluster mass function 
in Paper I, and the best estimate of Elmegreen et al. (2002) , 
2.2 < a < 2.4. 

Finally, in Fig. 5 we compare the clusters in the ring 
with those outside the ring. Within the observational un- 
certainties, we do not detect any differences in the slopes 
of either the age or the mass distributions. This is another 
argument highlighting the very young age of this cluster 
system compared to the expected characteristic disruption 
time-scale, even in the higher density environment of the 
starburst ring. 



3 STAR FORMATION IN THE BIRD'S HEAD 
GALAXY 

3.1 Observations and data preparation 

Optical observations of NGC 6745 in the minimum of four 
passbands spanning the entire wavelength range from the U 
(F336W) to the I band (F814W) were obtained between 18 
and 20 March 1996, as part of HST GO programme 6276 
(PI Westphal), using the Wide Field Planetary Camera 2 
(WFPC2). Multiple exposures were obtained through each 
filter to facilitate the removal of cosmic rays; the galactic 
centre was placed on the WF3 chip (pixel size 0.0997 arcsec) 
in all cases. An overview of the available HST observations 
is given in Table 1. 

Pipeline image reduction and calibration of the WFPC2 
images were done with standard procedures provided as part 
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Figure 4. (a) - The cluster formation rate (in number of clusters per Myr) in NGC 3310 as a function of age. The long-dashed line for 
log(Agc/yr) < 8 is the expected effect of a fading SSP, for a constant ongoing cluster formation rate, shifted vertically to best match the 
data points for log(Agc/yr) < 7.5. The short-dashed and dotted lines are the best-fitting and maximum deviating power-law slopes, as 
described in the text, (b) - Mass spectrum of the NGC 3310 clusters (number of clusters per mass bin); line encoding as in panel (a). 
The vertical dashed line at log(M c i/M0) = 4.4 corresponds to our completeness limit, for a cluster age of log(Age/yr) = 7.5, while the 
vertical dotted line corresponds to the same completeness limit at an age of log(Age/yr) = 8.0. 
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Figure 5. Comparison of the cluster formation rates and mass spectra of the clusters in vs. outside the NGC 3310 starburst ring. The 
age ranges in which stellar evolutionary fading and cluster disruption dominate are indicated; the vertical dash-dotted line in panel (b) 
indicates our completeness limit for a cluster age of log(Age/yr) = 7.5, while the vertical dotted line corresponds to the same completeness 
limit, but at an age of log(Age/yr) = 8.0. 



of the iraf/stsdas^ package, using the updated and cor- 



t The Image Reduction and Analysis Facility (IRAF) is dis- 
tributed by the National Optical Astronomy Observatories, which 
is operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National 



Science Foundation. STSDAS, the Space Telescope Science Data 
Analysis System, contains tasks complementary to the existing 
IRAF tasks. We used Version 2.2 (August 2000) for the data re- 
duction performed in this paper. 
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Table 1. Overview of the HST observations of NGC 6745 



Filter 


Exposure time 


ORIENT" 




(sec) 


(°) 


F336W 


2400, 7x2800 


145.52 


F555W 


2x1100 


145.52 




2x1300 


145.52 


F675W 


4x1300 


145.52 


F814W 


4x1300 


145.52 



NOTE: a — Orientation of the images (taken from the image header), mea- 
sured North through East with respect to the V3 axis (i.e., the X— Y diag- 
onal of the WF3 CCD +180°). 
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Figure 6. Formal completeness curves for the full WFPC2 field 
of view of the NGC 6745 observations. The different line styles 
refer to different passbands, as indicated in the figure. 

rected on-orbit flat fields and related reference files most 
appropriate for the observations. 

We obtained source lists and, subsequently, source pho- 
tometry for all of our NGC 6745 observations following pro- 
cedures identical to the ones described in detail in Paper I. 
For NGC 6745, at a distance D ~ 68 Mpc (m-M = 34.17; 
assuming a systemic velocity v r ,Q = 4545 ±60 km s _1 [Falco 
et al. 1999], a 208 km s _1 correction for infall of the Lo- 
cal Group towards the Virgo cluster, and Ho = 70 km s _1 
Mpc -1 ), all compact clusters and most star forming regions 
appear as point-like sources. The final list of verified cluster 
candidates detected at at least four times the background 
r.m.s. noise level in all four passbands contains 177 objects. 

We estimated the completeness of our source lists by 
using synthetic source fields consisting of PSFs. We created 
artificial source fields for input magnitudes between 20.0 and 
28.0 mag, in steps of 0.5 mag, independently for each of the 
F336W, F555W, and F814W passbands. We then applied 
the same source detection routines used for our science im- 
ages to the fields containing the combined galaxy image and 
the artificial sources. The results of this exercise are shown 
in Fig. 6. These formal completeness curves were corrected 
for the effects of blending or superposition of multiple ran- 
domly placed artificial PSFs as well as for the superposition 
of artificial PSFs on genuine objects. A detailed description 
of the procedures employed to obtain these completeness 
curves was given in Paper I. 

We found that the effects of image crowding are small: 
only < 1.5 — 2.5% of the simulated objects were not retrieved 
due to crowding, either in the artificial or in the combined 
frames. However, the effects of the bright and irregular back- 
ground and dust lanes are large, resulting in variable com- 



pleteness fractions across the galaxy images. As a general 
rule, however, the curves in Fig. 6 show that the formal 
completeness drops below ~ 90% for F336W > 23.5 mag, 
for F555W > 25.5 mag, and for F814W > 26 mag; the actual 
completeness limits are likely somewhat shallower, resulting 
from our selection criteria (cf. N. Bastian et al., in prep.). 

Foreground stars are not a source of confusion. The 
standard Milky Way star count models (e.g., Ratnatunga 
& Bahcall 1985) predict roughly 1-2 foreground stars for 
the equivalent standard filter of F555W < 24 in our field 
of view. Background objects may pose a (small) problem, 
however, in particular among the fainter sources, since we 
did not impose any roundness or sharpness constraints on 
our extended source detections, in order not to omit un- 
relaxed young clusters and star-forming regions from our 
final sample. However, such objects should be easily iden- 
tifiable once we have obtained aperture photometry for our 
complete source lists, as they are expected to have signif- 
icantly different colours. Background galaxies at redshifts 
greater than about 0.1 are expected to have extremely red 
(mF336W — wiF8i4w) colours compared to their local counter- 
parts and the young star clusters and star forming regions 
expected in NGC 6745. 

We will discuss the derived parameters and their impli- 
cations for the evolution of the galaxy's star cluster system 
in Sects. 3.3 and 3.4, respectively. 

Additionally, Hi 21cm observations (at ~ 1.399 GHz) 
were obtained with all 27 antennae of the Very Large Ar- 
ray (VLA) in C configuration, on 28 May 2000, with an 
on-source exposure time of 24765 s (~ 7 hr). Primary cali- 
brators were observed at the beginning and end; secondary 
(phase) calibrator integrations were uniformly interspersed 
throughout the run. The (u,v) plane coverage was almost 
perfectly isotropic, and 8192 CLEAN cycles (with uniform 
(u,v) weighting) yielded an essentially symmetric ~ 13.75 
arcsec FWHM beam. The 31 ~ 42.5 km s _1 channels were 
centered on 4541 km s _1 (channel 16), with channels 11 
through 21 containing all of the Hi emission. Continuum 
subtraction was accomplished in (u,v) space. We will dis- 
cuss these observations in detail in Sect. 3.2. 

3.2 The interaction geometry and its 
consequences 

The optical appearance of the NGC 6745 system, and in 
particular the locations of the numerous bright blue star 
forming complexes, are suggestive of a tidal passage by the 
small northern companion galaxy (NGC 6745c; nomencla- 
ture from Karachentsev, Karachentseva & Shcherbanovskii 
1978) across the eastern edge of the main galaxy, NGC 
6745a. The high relative velocities of the two colliding galax- 
ies likely caused ram pressures at the surface of contact be- 
tween the interacting interstellar clouds in both galaxies, 
which - in turn - are responsible for the triggering of en- 
hanced star formation, most notably in the interaction zone 
in between the two galaxies, NGC 6745b. 

We analysed our VLA Hi observations with as principal 
aim to assess whether we could confirm this scenario, at least 
qualitatively, and if we could, therefore, derive the interac- 
tion time-scale from our cluster analysis. The interaction- 
induced cluster formation scenario is at least qualitatively 
supported by the Hi zeroth-moment map of the system, 
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Figure 7. NGC 6745 Hi distribution shown as contours superimposed on a negative logarithmic reproduction of the F555W HST/WFPC2 
image of the galaxy. The VLA/C beam size is 13.75 arcsec (FWHM) and is indicated in the lower left corner. Contour intervals are 10 
per cent of the maximum flux density of the southern maximum, 9.565 X 10 2 Jy beam -1 m s -1 , except for the lowest contour, which is 
at 5 per cent of the maximum. Note the East- West asymmetry, the stretching of the lowest contours to the South, and the absence of 
detectable emission from the northern galaxy. For reference, the (J2000) centre coordinates of NGC 6745a and c are R.A. = 19 h 01 m 41.8 s , 
Dec = 40°44'40" and R.A. = 19 h 01 m 42.0 s , Dec = 40°45'35", respectively. 



shown overlaid on the optical WFPC2 image in Fig. 7, and 
the Hi velocity (first moment) map of Fig. 8. First, we point 
out that there is no evidence for Hi emission from the north- 
ern galaxy. This is consistent with the absence of significant 
recent star formation observed in this component (but see 
Sect. 3.3.2). 

The slight distortions in the symmetry of the overall 
Hi distribution along and perpendicular to the direction of 
elongation may, in fact, be consistent with tidal effects by a 
small, low- mass intruder (i.e., NGC 6745c). At the northern 
end of the Hi distribution the relevant time-scales are likely 
so short that we would not expect to see any significant 
amounts of tidally stretched Hi gas, which seems to be the 
case. At the southern end of the system, on the other hand, 
where the time elapsed since the alleged South-to-North pas- 
sage of the intruder is much longer, we see the expected 
(asymmetric) extension of the lower contours of the Hi dis- 
tribution, as well as a high velocity of approach at the very 
bottom of our field of view, which may be qualitatively con- 
sistent with a trajectory of approach for the intruder. This 
interpretation is quantitatively supported by CO(l — 0) line 
observations of this system by Zhu et al. (1999), who find 
a similar velocity gradient between the southern component 
and the centre of the main galaxy. 

The velocity data (Fig. 8) were fit with a Brandt func- 



tion (Brandt 1960, Brandt & Belton 1962), weighting the 
velocities with the square of the intensity data (Fig. 7). The 
residual map showed excellent agreement over all but the 
lowest-intensity regions. The best-fitting parameters include 
(i) systemic velocity, v r .Q = 4528 km s~ at 4 arcsec West 
and 2 arcsec South of the optical nucleus; (ii) position an- 
gle of the line of nodes = 36 degrees; (iii) inclination of the 
fundamental plane of the galaxy i = 47 degrees; and (iv) 
half- amplitude value = 185 km s~ reached at ~ 32 arc- 
sec from the location of the kinematic centre. The sense of 
the rotation of the Hi means that the perturbing galaxy, 
NGC 6745c, is in a prograde orbit, which favours a strong 
and rapid gravity-wave response in the affected galaxy, NGC 
6745a. 



VLA continuum observations at 1.425 GHz (Condon et 
al. 1996) show a more compact morphology, which more 
closely follows the system's stellar distribution and far- 
infrared appearance (IRAS; e.g., Bushouse et al. 1988, Con- 
don et al. 1995, Sanders et al. 1995). As opposed to our 
~ 1.4 GHz line observations that predominantly trace the 
cold Hi component, the continuum flux has a non-thermal 
origin, as evidenced by the 1.4 — 5 GHz spectral index of 
between +0.75 and +0.83 (Condon et al. 1991, 1995). 
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Figure 8. NGC 6745 Hi radial velocity distribution shown as contours superimposed on a negative logarithmic reproduction of the 
F555W HST/WFPC2 image of the galaxy. The VLA/C beam size is 13.75 arcscc (FWHM) and is indicated in the lower left corner. 
Contours are at intervals of 20 km s _1 beginning at 4400 km s _1 for the smallest entire contour South West of the center enclosing the 
minimum of 4388 km s~ 1 and ending with the smallest entire contour North East of the center enclosing the maximum of 4661 km s _1 . 
The "boxy" outer line encloses map data > 8.14 X 10 -4 Jy beam -1 (3 a) used in the moment calculations. The contours exhibit, in 
general, the splayed structure that is characteristic of rotating disks common in late-type spiral galaxies. 



3.3 Properties of the star- forming regions 

Since the morphology and appearance of the NGC 6745 sys- 
tem are consistent with the scenario favouring a very recent 
tidal encounter and its associated triggered star and star 
cluster formation, we will now use the main properties of the 
star clusters and star forming complexes, such as their ages, 
masses and metallicities, to derive boundary conditions for 
the physical conditions governing such minor merger events. 

With the lessons learnt from Paper I in mind, we applied 
a similar three-dimensional \ 2 minimisation (with respect to 
the Anders & Fritze-v. Alvensleben [2003] models) to the 
spectral energy distributions (SEDs) of our star cluster can- 
didates and star forming complexes to obtain the most likely 
combination of age, mass, metallicity Z and internal extinc- 
tion E(B — V) (assuming a Calzetti et al. [2000] starburst 
galaxy-type extinction law for the internal extinction) for 
each individual object. Galactic foreground extinction to- 
wards NGC 6745 was taken from Schlegel et al. (1998). The 
resulting age, mass, metallicity and extinction distributions 
for the NGC 6745 objects are shown in Fig. 9. 



Fig. 9a) might be an artefact caused by the age-metallicity 
degeneracy. Therefore, we display these particular clusters 
as the shaded histograms in Figs. 9b, c and d. It is clear from 
an examination of the NGC 6745 cluster metallicity distri- 
bution in Fig. 9c, that the clusters in the young-age peak 
indeed compose a higher-than-average metallicity peak, thus 
highlighting the caution required to avoid confusion arising 
from this degeneracy. As a consequence of the underesti- 
mated ages for the clusters affected by the age-metallicity 
degeneracy, their masses are also biased towards lower-than- 
average masses. If we now turn this argument around, and 
examine the properties of the clusters in the second metal- 
licity peak at 0.015 < Z < 0.020, we find that those clus- 
ters do indeed dominate the young-age peak (as expected; 
not shown), and are generally biased towards younger ages 
overall. This confirms the importance of the age-metallicity 
degeneracy for the broad-band filters used in this paper, and 
for the young ages of the star clusters and star-forming com- 
plexes in NGC 6745. The mass and extinction properties of 
the clusters in this secondary metallicity peak are consistent 
with those of the full cluster sample. 



3. 3. 1 Effects of the age-metallicity degeneracy 

Based on the knowledge gained in Paper I, we suspected 
that the strong peak at young ages (log(Age/yr) ~ 6.8; see 
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Figure 9. Age, mass, metallicity and extinction distributions of the clusters in NGC 6745, based on our 4-passband fitting technique used 
to determine age, metallicity and extinction estimates simultaneously. Although the main properties of these distributions arc reasonably 
correct, the shaded histograms are affected significantly by the age— metallicity degeneracy: the shaded histograms, corresponding to the 
young-age peak in panel (a), show a distinct second peak at close-to-solar metallicity (panel c), compared to the other clusters that peak 
at much lower metallicities. The shaded histogram in the inset in panel (c) is the metallicity distribution of the clusters in the most 
actively star forming region connecting the main galaxy with its smaller companion. 

fitting routines to the observational data, assuming a more 
appropriate fixed metallicity of Z — 0.008 for all clusters. 
The resulting age and mass distributions, are shown in Figs. 
10b and c. We also show (in Fig. 10a) the clusters' locations 
in the (age vs. mass) plane with respect to our ~ 90% com- 
pleteness limit at w-fsssw — 23.5 mag (imposed by the shal- 
low depth of the F336W observations; cf. Fig. 6). The new 
distribution of extinction estimates, based on the Z — 0.008 
assumption, does not differ significantly from that of Fig. 
9d. 

[We note that, while this procedure likely results in the 
most realistic cluster parameters for the NGC 6745 system 
(due to the availability of only four optical passbands), the 
situation for NGC 3310 is different. Because of the availabil- 
ity of observations of the NGC 3310 cluster system cover- 
ing a larger number of passbands and a greater wavelength 



3.3.2 Case for a fixed mean metallicity 

In the inset of Fig. 9c we show the metallicity distribution of 
the entire cluster sample once again, with the clusters from 
the interaction zone in between the two galaxies highlighted 
as the inner, shaded histogram. This shaded histogram is 
characterized by two clear peaks in metallicity. However, if 
these clusters were all formed approximately coevally (as 
implied by their derived ages) and in a tightly confined spa- 
tial area, their metallicities must be very similar; a clearly 
double-peaked distribution is therefore physically unrealis- 
tic and therefore most likely another artefact caused by the 
age-met allicity degeneracy. 

Based on these arguments, it is most likely that the 
NGC 6745 young cluster system is best approximated by a 
metallicity consistent with that of the low-metallicity peak 
in Fig. 9c, at Z < 0.015. Therefore, we reapplied our SED 
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Figure 10. (a) - Distribution of the NGC 6745 clusters in the (age vs. mass) plane, assuming a fixed metallicity of Z = 0.008 for all 
clusters. The solid bullets represent the clusters in the most actively star forming region, i.e. NGC 6745b. Overplotted is the expected, 
age-dependent detection limit, predicted by the Anders & Fritze-v. Alvensleben (2003) models for the appropriate stellar IMF. This 
model prediction is based on a detection limit of mp33gw = 23.5 and (m — M)jmqci6745 = 34.17, assuming no extinction, (b) and (c) - 
Distributions of, respectively, the ages and masses of the compact clusters in NGC 6745. The shaded histograms show the distributions 
of the clusters in NGC 6745b. The uncertainties are of the order of the histogram bin sizes. 



range, we concluded that fitting all of our free parameters 
simultaneously resulted in the best parameter estimates.] 

There is an enhancement of cluster formation around 
log(Age/yr) ~ 7, which is therefore our best estimate 
of the interaction time-scale. The clusters younger than 
log(Age/yr) ~ 7.4 (~ 25 Myr) are all located in either the 
interaction zone NGC 6745b (see the shaded histogram in 
Fig. 10b), or among the sprinkling of blue objects resembling 
a spiral arm on the eastern edge of the main galactic disc. 
This is consistent with tidally-triggered star (cluster) forma- 
tion along the proposed trajectory of the companion galaxy 
across the main disc component of NGC 6745 (see Sect. 3.2). 
The lower age cut-off at log(Age/yr) = 6.6 (4 Myr) is ar- 
tificial, and corresponds to the youngest ages included in 
our SSP models (Schulz et al. 2002, Anders & Fritze-v. Al- 
vensleben 2003) ; a fraction of our NGC 6745 cluster sample 



may therefore have even younger ages. Thus, it seems likely 
that the interaction began to have a significant effect on the 
star (cluster) formation history some 25 Myr ago, and has 
been ongoing until the present. 

We detected three compact clusters that are clearly as- 
sociated with the companion galaxy itself; they are unlikely 
to originate from the interaction-induced starburst region, 
NGC 6745b, based on their large projected distances from 
the interaction zone. Two of these have ages of ~ 100 and 
~ 140 Myr, and are located in the general field of its disc. 
The third, located near the centre of the companion galaxy, 
is of similar age as the starburst in NGC 6745a and b, ~ 8 
Myr. Close inspection of the overall colour distribution of 
the field star population in NGC 6745c reveals a bluer cen- 
tre, which might therefore imply that the tidal interaction 
with the main NGC 6745 galaxy has also induced enhanced 
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star formation in the companion galaxy, but at a much lower 
level due to the absence of significant amounts of gas (see 
Sect. 3.2). The little available gas in NGC 6745c might not 
have been sufficient for the formation of star clusters. 

There is no clear spatial dependence of the cluster 
masses, although the lowest-mass clusters appear to origi- 
nate in the lowest-density regions. This is most likely an ob- 
servational selection effect, however. Although the galaxy's 
distance, combined with the photometric depth of the obser- 
vations, prevents us from probing below log(M c i/M ) = 5.0 
for clusters older than ~ 10 Myr, we point out that there is a 
significant population of high-mass clusters or cluster com- 
plexes, with masses in the range 6.5 < log(M c i/M ) < 8.0. 
We have checked the reality and signal-to-noise ratio of the 
SED of each of these objects individually, and have only re- 
tained those objects that appear real beyond any reasonable 
doubt. These clusters do not have counterparts among the 
Galactic GCs (e.g., Mandushev, Staneva & Spassova 1991, 
Pryor & Meylan 1993). In fact, these masses are similar to 
or exceed the spectroscopically confirmed mass estimates of 
the so-called "super star clusters" (SSCs) in M82 (M82-F; 
Smith & Gallagher 2001), and the Antennae galaxies (Men- 
gel et al. 2002). Our detection of similarly massive SSCs in 
NGC 6745, which are mostly located in the intense inter- 
action zone, NGC 6745b, supports the scenario that such 
objects form preferentially in the extreme environments of 
interacting galaxies. 

We caution, however, that these massive SSC candi- 
dates may not be gravitationally bound objects, but more 
diffuse star forming regions or aggregates of multiple unre- 
solved clusters instead. Even with the unprecedented HST 
spatial resolution, we cannot distinguish between these pos- 
sibilities. Nevertheless, we measure an effective radius for the 
most massive object (Ma ~ 5.9 x 1O 8 M ) of only 7? ef f ~ 16 
pc. However, this object appears very elongated, or may in 
fact be a double cluster. We should keep in mind, of course, 
that this high mass estimate is a strong function of the (low) 
metallicity assumed; if we had assumed solar metallicity for 
this object, the derived age would have been significantly 
smaller (~ 10 — 20 Myr vs. ~ 1 Gyr), and the mass could 
be smaller by a factor of > 10. Even so, if we could confirm 
this mass estimate spectroscopically, either of the subcom- 
ponents would be the most massive cluster known to date, 
significantly exceeding cluster W3 in NGC 7252, which has a 
mass of about (3 — 18) xlO 7 M , depending on the age, metal- 
licity and IMF assumed (Schweizer & Seitzer 1998, Maraston 
et al. 2001). 

3.4 The initial cluster mass distribution 

As for NGC 3310, in Figs. 11a and b we show the NGC 6745 
cluster formation rate (per linear time period) and its mass 
spectrum, respectively. For our analysis of the cluster disrup- 
tion time-scales in NGC 6745 we will make the same assump- 
tions as described in Sect. 2.2. The long-dashed line in Fig. 
11a for the youngest ages, log(Age/yr) < 8, is the expected 
effect of a fading SSP in the F336W filter (which determines 
our sample completeness for this galaxy) , for a constant on- 
going cluster formation rate, shifted vertically to best match 
the data points for log(Age/yr) < 7.25, and limited by the 
completeness in the F336W observations. For this band, the 
predicted fading line is horizontal at log(Age/yr) < 6.5 and 



has a slope of —1.12 for older clusters, in the parameter space 
of Fig. 11a. The short-dashed lines in both Figs. 11a and b 
are the best-fitting power-law slopes, obtained for the age 
and mass ranges log(Age/yr) > 7.4 and log(M c i/M ) > 5.6, 
while the dotted lines indicate the likely uncertainties in the 
slopes. The latter were obtained from subsets of the data 
points used for the dashed fits, which were determined over 
the full ranges. The vertical dashed and dotted lines in Fig. 
lib at log(M c i/M ) = 4.75 and log(M c i/M ) = 5.4 corre- 
spond to our formal completeness limit, for a peak young 
cluster age of log(Age/yr) = 7.0 and for the significant frac- 
tion of sample clusters with ages up to log(Age/yr) =8.0, 
respectively. Thus, it follows that incompleteness plays at 
least some role for masses up to log(M c i/M ) ~ 5.4. 

For log(Age/yr) > 7.25 we find a best-fitting slope of 
(1 - a)/7 = -1.86(±0.20)t2 og, where the first uncertainty 
represents the formal uncertainty in the fit and the second 
the likely uncertainty resulting from our choice of fitting 
range (i.e., the dotted fits in Fig. 11a). Assuming that a = 
2.0 (see below) we find that 7 = 0.54 ± 0.06, close to the 
mean value of 0.62 derived by BL03 for cluster samples in 
four galaxies. 

The crossing point of the two power laws of the age dis- 
tribution is at log(t croS s/yr) — 7.3. Substituting these values 
into Eq. (15) of BL03, we obtain a characteristic disrup- 
tion time-scale for 1O 4 M clusters of log(4 is /yr) = 7 - 75 - 
This implies that the disruption time-scale of clusters more 
massive than ~ 4 x 1O 5 M0 is longer than 5.5 x 10 8 yr. 
All clusters in our sample, with one exception, are younger 
than this. The mass distribution of clusters more mas- 
sive than 4 x 1O 5 M , shown in Fig. 11a, has a slope of 
— 1.96(±0.15) ± 0.19, where the first uncertainty again rep- 
resents the formal uncertainty in the fit and the second the 
likely uncertainty resulting from our choice of fitting range 
(i.e., the dotted fits in Fig. lib). Since these massive clus- 
ters have not yet suffered disruption, the measured slope 
represents that of the cluster IMF of NGC 6745. 



4 THE INITIAL CLUSTER MASS FUNCTION 

For star cluster systems as young as the tidally-induced clus- 
ter populations in NGC 3310 and NGC 6745, for which the 
age of the starburst in which they were formed is signif- 
icantly smaller than their respective characteristic cluster 
disruption time-scales, we have shown that the application 
of the empirical cluster disruption models of BL03 results in 
an independent estimate of the initial cluster mass function 
slope, a. For the NGC 3310 and NGC 6745 cluster systems, 
a = 2.04(±0.23)±oi3 and 1.96(±0.15) ±0.19, respectively, 
for masses Ma > 1O 5 M and M cl >4x 1O 5 M . 

The slopes of the CLFs of young star cluster systems, 
and to some lesser extent those of their mass functions as 
well, have been studied extensively - in particular since the 
launch of the HST - in a wide variety of star-forming envi- 
ronments, including interacting and starburst galaxies, cir- 
cumnuclear starburst rings, and also in "normal" spiral and 
irregular galaxies by means of the Hn region luminosity func- 
tion (regarding the latter, see e.g. Kennicutt & Hodge 1980, 
Kennicutt et al. 1989, Arsenault et al. 1990, Cepa & Beck- 
man 1990, Caldwell et al. 1991, Banfi et al. 1993, Elmegreen 
& Salzer 1999; and simulations by Oey & Clarke 1998). In 
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Figure 11. (a) — The cluster formation rate (in number of clusters per Myr) in NGC 6745 as a function of age. The long-dashed line for 
log(Age/yr) < 8.0 is the expected effect of a fading SSP, for a constant ongoing cluster formation rate, shifted vertically to best match 
the data points for log(Age/yr) < 7.25, and limited by the completeness of our F336W observations. The short-dashed and dotted lines 
are our best fit for greater ages and the maximum likely deviation from this best fit, respectively, (b) - Mass spectrum of the NGC 6745 
clusters (number of clusters per mass bin); line encoding as in panel (a). The vertical dashed line at log(Af c i/ATQ) = 4.75 corresponds to 
our 90% completeness limit for a cluster age of log(Agc/yr) = 7.0. The vertical dotted line corresponds to the same completeness limit, 
but for a cluster age of log(Age/yr) = 8.0. 



Table 2, we give an overview of the CLF (and, in a few 
cases, of the mass function) slopes based on determinations 
available in the literature, both for young cluster systems in 
interacting and starburst galaxies, and for cluster popula- 
tions in circumnuclear starburst rings. 

In addition to these galaxies, Meurer et al. (1995) found 
that the CLF slope of the combined sample of SSCs in 9 
starburst galaxies, for Mf220W J; — 14, is consistent with 
a — 2; Maoz et al. (1996) constructed the combined CLF 
for the circumnuclear ring clusters in NGC 1433, NGC 1512 
and NGC 2997, based on UV observations, and reached a 
similar conclusion for the high-luminosity end of their CLF. 
Finally, Crocker et al. (1996) analysed the Hll region lumi- 
nosity function in 32 (pseudo-)ringed SO - Sc galaxies, and 
found a mean (a) = 1.95 ± 0.25. 

Thus, in all but a few cases (e.g., M82 B), the CLF 
and, where derived, the mass function slopes of the young 
clusters are similar (a ~ 2) among a wide variety of galax- 
ies, including those in which the star and cluster formation 
process is dominated by strong tidal forces, and systems of 
a more quiescent nature (e.g., normal spiral and irregular 
galaxies; Kennicutt et al. 1989, Banfi et al. 1993, Elmegreen 
& Salzer 1999) . One should keep in mind that the observed 
CLF is the result of both the original mass distribution of the 
young star clusters and the giant molecular clouds (GMCs) 
they originated from, and subsequent cluster disruption pro- 
cesses (e.g., Harris & Pudritz 1994, Bik et al. 2003, BL03, 
de Grijs et al. 2003b); in de Grijs et al. (2003b) and Par- 
mentier et al. (2003) we showed that the M82 B cluster 
system must have undergone significant disruption owing to 
the unusually short characteristic cluster disruption time- 
scale, log(t dis /yr) ~ 7.5 + 0.62 x log(M c i/lO 4 M ) (de Grijs 



et al. 2003b) , so that the observed shallow(er) CLF and mass 
function slopes do not represent the initial cluster mass func- 
tion slope anymore. However, in most other young cluster 
systems known the expected characteristic cluster disrup- 
tion time-scales are well in excess of the median age of the 
cluster system, so that at least the high-luminosity (high- 
mass) end of the CLF (and cluster mass function) remains 
likely unaffected by significant alterations caused by disrup- 
tion (e.g., Buta et al. 1999, Zhang & Fall 1999, Bik et al. 
2003, Paper I). In this paper, we have derived very similar 
mass function slopes, and have also properly accounted for 
the effects of cluster disruption. 

Under the usual assumption that the stellar mass-to- 
light ratio does not vary significantly over the age range of a 
given young cluster system, or - alternatively - after correct- 
ing the observed CLF to a common age (e.g., Meurer 1995, 
Fritze-v. Alvensleben 1999, de Grijs et al. 2001, 2003a,b), 
the resulting CLF slope can be interpreted in terms of the in- 
trinsic cluster mass distribution. The observed mass function 
slope, a ~ 2, is very similar to those of the mass functions of 
Galactic open star clusters (e.g., van den Bergh & Lafontaine 
1984, Elson & Fall 1985), GMCs (e.g., Casoli, Combes & 
Gerin 1984, Sanders, Scoville & Solomon 1985, Solomon et 
al. 1987, Solomon & Rivolo 1989, Williams & McKee 1997) 
and GMC cores (e.g., Harris & Pudritz 1994, Williams, de 
Geus & Blitz 1994, McLaughlin & Pudritz 1996 and refer- 
ences therein; see also the review by Blitz 1991), and to that 
of the young compact cluster system in the Large Magellanic 
Cloud (LMC; e.g., Elson & Fall 1985, Elmegreen & Efremov 
1997). Thus, it appears that the formation of young compact 
star clusters, such as the ones discussed in the current series 
of papers, and of open clusters, GMCs and stellar OB asso- 
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Table 2. Cluster luminosity and mass function slopes based on previous determinations in the literature 



Galaxy 



Slope ( — a) 



Reference 



NGC 1275 

NGC 1316 

NGC 3256 
NGC 3610 

NGC 3921 

NGC 4038/39 



NGC 7252 

M51 
M82 B 



NGC 1097 
NGC 1512 
NGC 2997 

NGC 3310 



NGC 5248 
NGC 6951 
ESQ 565-11 



~ -2 

f 2 

-1.7 ± 0.1 

~ -1.8 
-1.78 ± 0.05 
-1.90 ± 0.07 
-2.1 ± 0.3 
-2.01 ± 0.22 
-1.78 ± 0.05 
-2.12 ± 0.04 
-2.6 ± 0.2 
-1.7 ± 0.2 
-1.95 ± 0.03 
-2.00 ± 0.08 

~ -2 
-1.90 ± 0.04 
-1.80 ± 0.07 
-2.1 ± 0.3 
-2.0 ± 0.05 
-1.2 ± 0.3 
-1.4 ± 0.2 



r 2 

~ -2 

2 

-1.97 ± 0.03 
-2.2 ± 0.03 
-2.4 ± 0.04 
-1.75 ± 0.03 
-1.8 ± 0.4 

2 

-2.12 ± 0.04 
-2.18 ± 0.06 



A. Interacting and Starburst Galaxies 
My < -12 

M B < -8 

3 Gyr-old merger remnant 

M B < -9.5; Mj < -10.5 
V < 26; red, metal-rich clusters 
corrected for observational scatter 
My < 8.5 

corrected for contamination by old GCs 

— 15.5 < My < —10; no completeness corrections 
— 14 < My < —8, but contaminated by stars 

— 14 < My < —10.4, corrected for contamination by stars 

— 10.4 < My < —8, corrected for contamination by stars 
6.4 < log(Age/yr) < 6.8; 10 4 < M cl < 1O 6 M 

7.4 < log(Age/yr) < 8.2; 10 4 < M cl < 1O 6 M 
My < -12 

My < —8, inner sample 
My < —8, outer sample 

2.5 X 10 3 < M cl < 5 X 1O 4 M ; log(Age/yr) < 7 

2.5 X 10 3 < M cl < 2 X 1O 4 M ; log(Age/yr) < 7 
My < -11 

8.6 < log(Age/yr) < 9.1 

B. Circumnuclcar Starburst Rings 
-14 < My < -11 



10 3 < M cl < 10 5 M, 







My < -9 

M B < -10.5 

from HST X-band data 

large-scale cluster complexes, K band 

17.7 < m F606w < 20.2 

10 3 < M cl < 1O 5 M 

My < -9 

My < -8.9; Cluster ages 4-6 Myr 



Holtzman et al. (1992) 
Carlson et al. (1998) 
Goudfrooij et al. (2001) 
Whitmorc et al. (2002) 
Zepf et al. (1999) 
Whitmorc et al. (2002) 

Schweizer et al. (1996) 

Whitmore & Schweizer (1995) 
Whitmorc et al. (1999) 



Zhang & Fall (1999) 

Whitmorc et al. (1993) 
Miller et al. (1997) 

Bik et al. (2003) 

do Grijs et al. (2001) 
Parmcnticr et al. (200.3) 



Barth et al. (1995) 
Maoz et al. (2001) 
Maoz et al. (1996) 
Elmegrecn et al. (1999) 
Elmegreen et al. (2002) 



Paper I 

Maoz et al. (2001) 
Elmegrecn et al. (1999) 
Buta et al. (1999) 



ciations (e.g., Elmegreen & Efremov 1997, Elmegreen 2002) 
is driven by similar physical processes, which results natu- 
rally from models advocating that the mass distribution of 
the progenitor GMCs is well approximated by a turbulence- 
driven fractal structure, allowing cloud growth by agglom- 
eration (e.g., Harris & Pudritz 1994, Elmegreen & Efremov 
1997, Elmegreen et al. 1999, Elmegreen 2002). 



5 SUMMARY AND CONCLUSIONS 

Both NGC 3310 and NGC 6745 are examples of nearby in- 
teracting galaxies undergoing intense, tidally-triggered star- 
bursts. The production of luminous, compact star clusters 
seems to be a hallmark of intense star formation. It is likely 
that a large fraction of the star formation in starbursts takes 
place in the form of such concentrated clusters. The basic 
cluster properties, including their sizes, luminosities, and 
- in several cases - spectroscopic mass estimates are en- 
tirely consistent with what is expected for young Galaxy- 
type GCs. Young compact star clusters are therefore impor- 
tant because of what they can tell us about GC formation 
and evolution. They are also important as probes of the his- 
tory of star formation, chemical evolution, IMF, and other 
physical characteristics in starbursts. 

In this paper we have used the basic properties of the 
rich young star cluster systems in both starburst galaxies, 
including estimates of the individual cluster ages, masses, 
metallicities, and extinction values, to derive the cluster for- 



mation history and the subsequent evolution of the star clus- 
ter systems as modified by the effects of cluster disruption. 
Cluster disruption processes must be taken into account for 
the determination of the cluster formation history from the 
age distribution of a magnitude-limited cluster sample, be- 
cause the observed age distribution is that of the surviving 
clusters only. 

For the NGC 3310 star cluster system we used the ba- 
sic properties obtained in Paper I, which were based on the 
detailed comparison of updated model SSPs (including the 
contributions of gaseous, nebular emission) with the clusters' 
SEDs derived from broad-band HST imaging observations. 
We applied our improved understanding of the systematic 
uncertainties involved in such an exercise (Paper I) to simul- 
taneously determine robust cluster ages, masses, metallici- 
ties and extinction values for the tidally-induced young star 
cluster system in NGC 6745, based on imaging observations 
in the minimum of four broad-band passbands covering the 
entire optical wavelength range from the HST-equivalent U 
to /-band filters. We further expanded our analysis of the 
systematic uncertainties involved in this type of analysis by 
examining the effects of a priori assumptions of the individ- 
ual cluster metallicities. We find that - if observations cov- 
ering a sufficiently large optical/NIR wavelength range are 
available - (i) it is very important to determine all of the free 
parameters (age, metallicity, and extinction, and the corre- 
sponding mass estimates) independently for each individual 
cluster, instead of assuming a generic value for any of these 
parameters (such as assuming a fixed metallicity); and (ii) 
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the widths of the resulting age and mass distributions, un- 
der the assumption of a fixed metallicity for all clusters, are 
broader than their intrinsic widths owing to the propagation 
of model and measurement uncertainties. 

The central morphology of NGC 3310 is dominated by 
a bright starburst ring containing a large number of young 
star clusters. The age (and metallicity) distributions of both 
the clusters in the ring and of the sample of clusters out- 
side the ring are statistically indistinguishable, with a peak 
in the age distribution at log(Age/yr) ~ 7.45, but the ring 
clusters appear to peak at slightly greater masses than those 
found outside the ring. The overabundance of low-mass clus- 
ters outside the ring compared to their counterparts in the 
starburst ring might (at least partially) be due to a spa- 
tial dependence of the completeness fraction, but there is 
a clear and significant excess of higher-mass clusters in the 
ring compared to the non-ring cluster sample, even after tak- 
ing the systematic uncertainties in our mass estimates into 
account. It appears that the physical conditions in the star- 
burst ring, such as caused by the effects of the proposed bar- 
driven instabilities, the higher density of the ISM, and the 
associated higher likelihood to encounter significant propa- 
gating shock waves, may be conducive for the formation of 
higher-mass star clusters, on average, than in the relatively 
more quiescent environment of the main galactic disc. 

We furthermore conclude that the starburst producing 
the NGC 3310 ring clusters has been ongoing for at least the 
past 40 Myr, at an approximately constant cluster formation 
rate during the burst period. Star cluster formation has pro- 
ceeded at a similar rate in the central disc of the galaxy; due 
to the extreme youth of both samples of starburst-induced 
star clusters, statistically significant differences between the 
age distributions of the ring and non-ring populations have 
not yet had time to develop, which implies that we are still 
observing at least part of the cluster sample in the environ- 
ment determined by their initial conditions. 

For the NGC 6745 cluster system we derive an median 
age of ~ 10 Myr, assuming a fixed metallicity of Z — 0.008 
to avoid the significant effects of the age-met allicity degen- 
eracy for these young ages. Based on the age distribution of 
the star clusters, and the Hi morphology of the interacting 
system as observed with the VLA, we qualitatively confirm 
the interaction-induced enhanced star formation scenario. 
NGC 6745 contains a significant population of high-mass 
clusters, with masses in the range 6.5 < log(M c i/M Q ) < 8.0. 
These clusters do not have counterparts among the Galac- 
tic GCs, but are similar to or exceed the spectroscopically 
confirmed mass estimates of the SSCs in M82 and the Anten- 
nae galaxies. Our detection of such massive SSCs in NGC 
6745, which are mostly located in the intense interaction 
zone, NGC 6745b, supports the scenario that such objects 
form preferentially in the extreme environments of interact- 
ing galaxies. 

For star cluster systems as young as the tidally-induced 
cluster populations in NGC 3310 and NGC 6745, for which 
the age of the starburst in which they were formed is sig- 
nificantly lower than their respective characteristic cluster 
disruption time-scales (of, respectively, log(£4 ls /yr) = 8.05 
and 7.75), we have shown that the application of the em- 
pirical cluster disruption models of BL03 results in an inde- 
pendent estimate of the initial cluster mass function slope, 
a. For the NGC 3310 and NGC 6745 cluster systems, 



a = 2.04(±0.23)tg;Jij and 1.96(±0.15) ±0.19, respectively, 
for masses M c \ > 10 5 M Q and M c i > 4 x 10 5 M Q . These mass 
function slopes are consistent with those of other young star 
cluster systems in interacting and starburst galaxies, cir- 
cumnuclear starburst rings, and of the Hn region luminosity 
(and mass) functions in "normal" spiral and irregular galax- 
ies. They are also very similar to the mass function slopes of 
Galactic open clusters and (OB) associations, of GMCs and 
their cores, and of the young compact clusters in the LMC. It 
has been shown previously that this can be understood if the 
mass function of the progenitor GMC is approximated by a 
turbulence-driven fractal structure, allowing cloud growth 
by agglomeration. 
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